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The critical behaviour of systems belonging to the three-dimensional Ising universality class is studied theo- 
retically using the collective variables (CV) method. The partition function of a one-component spin system is 
calculated by the integration over the layers of the CV phase space in the approximation of the non-Gaussian 
sextic distribution of order-parameterfluctuations(the p 6 model). A specific feature of the proposed calculation 
consists in making allowance for the dependence of the Fourier transform of the interaction potential on the 
wave vector. The inclusion of the correction for the potential averaging leads to a nonzero critical exponent of 
the correlation function r\ and the renormalization of the values of other critical exponents. The contributions 
from this correction to the recurrence relations for the p 6 model, fixed-point coordinates and elements of the 
renormalization-group linear transformation matrix are singled out. The expression for a small critical exponent 
n is obtained in a higher non-Gaussian approximation. 
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1. Method 

We shall use the approach of collective variables (CV) , which allows us to calculate the expres- 
sion for the partition function of a system and to obtain complete expressions for thermodynamic func- 
tions near the phase-transition temperature T c in addition to universal quantities (i.e., critical exponents). 
The CV approach is non-perturbative and similar to the Wilson non-perturbative renormalization-group 
(RG) approach (integration on fast modes and construction of an effective theory for slow modes) 1 3- 5] . 

The term collective variables is a common name for a special class of variables that are specific for 
each individual physical system. The CV set contains variables associated with order parameters. Conse- 
quently, the phase space of CV is most natural in describing a phase transition. For magnetic systems, the 
CV pk are the variables associated with the modes of spin-moment density oscillations, while the order 
parameter is related to the variable po, in which the subscript "0" corresponds to the peak of the Fourier 
transform of the interaction potential. The methods available at present, make it possible to calculate 
universal quantities to a quite high degree of accuracy (see, for example, @]). The advantage of the CV 
method lies in the possibility to obtain and analyse thermodynamic characteristics as functions of micro- 
scopic parameters of the original system. The use of the non-Gaussian basis distributions of fluctuations 
in calculating the partition function of a system does not bring about a problem of summing various 
classes of divergent (with respect to the Gaussian distribution) diagrams at the critical point. A considera- 
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tion of the increasing number of terms in the exponent of the non-Gaussian distribution is an alternative 
to the use of a higher-order perturbation theory based on the Gaussian distribution. 

The integration of partition function begins with the variables p^ having a large value of the wave 
vector k (of the order of the Brillouin half-zone boundary) and terminates at pk with k —> 0. For this 
purpose, we divide the phase space of the CV pk into layers with the division parameter s. In each nth 
layer (corresponding to the region of wave vectors B n+ \ < k s£ B n , B n+ \ - B n /s, s > 1), the Fourier trans- 
form of the interaction potential is replaced by its average value. This simplified procedure leads to a 
zero value of the critical exponent 77 characterizing the behaviour of the pair-correlation function at the 
critical temperature T c . 

2. The setup 

The object of investigation is a three-dimensional (3D) Ising-like system. The Ising model, despite its 
simplicity, has, on the one hand, a wide scope of realistic applications, and, on the other hand, it can be 
considered as a model, which serves as a standard in studying other models possessing a much more 
complicated construction. 

In our previous calculations (for example, I7H10I1). we assumed that the correction for the potential 
averaging is zero. As a result, we lost some information in the process of calculating the partition function 
of the system. In particular, the critical exponent r\ was equal to zero. 

In lllll , the correction for the averaging of the Fourier transform of the potential is taken into account 
in the simplest non-Gaussian approximation (the p 4 model based on the quartic fluctuation distribution). 
The inclusion of this correction gives rise to a nonzero value of the critical exponent 77. Recurrence rela- 
tions (RR) between the coefficients of the effective distributions take another form as compared with the 
case when 77 = 0. A fixed point is shifted. Critical exponents of the correlation length v, susceptibility 7 
and specific heat a are renormalized. Critical amplitudes are also modified. As is seen from table [1] the 
inclusion of a nonzero exponent 77 within the CV method reduces the critical exponent v (like in the non- 
Table 1. Estimates of the critical exponents for the p 4 model and the RG parameter s = 4 in the case of the 
correction for the potential averaging not taken into account (AO(A;) = 0) and in the case of the correction 
for the potential averaging taken into account (AO(fc) 4- 0). 



Condition 


n 


V 


r 


a 


A<P(fc) = 





0.612 


1.225 


0.163 


AO(fc) *0 


0.024 


0.577 


1.141 


0.268 



perturbative RG approach (ill)- m order to obtain better quantitative estimates of v and other critical 
exponents, it is necessary to use the distributions of fluctuations more complicated than the quartic dis- 
tribution. In the case of 77 = 0, the critical exponent of the correlation length for these distributions takes 
on larger values than v for the p 4 model (figure[l} |2, 8]. The results of calculations and their comparison 
with the other authors' data show that the sextic distribution (the p 6 model) provides a more adequate 
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Figure 1. Evolution of the critical exponent of the correlation length v with an increasing parameter of 
division of the CV phase space into layers s. Curves 1, 2, 3 and 4 correspond to the p 4 , p 6 , p 8 and p 10 
models, respectively. 
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quantitative description of the critical behaviour of a 3D Ising ferromagnet than the quartic distribu- 
tion [2j[l3j]. The sextic distribution for the modes of spin-moment density oscillations is presented as an 
exponential function of the CV whose argument includes the sixth power of the variable in addition to 
the second and the fourth powers. 

In the present publication, our aim is to investigate the effect of the above-mentioned correction for 
the potential averaging on the critical properties of a 3D Ising-like system and to elaborate a technique for 
calculating the small critical exponent rj in the p 6 model approximation. The analytic results, obtained in 
this higher non-Gaussian approximation, provide the basis for accurate analysis of the behaviour of the 
system near T c with allowance for the exponent r\. 

The developed approach permits to perform the calculations for a one-component spin system in real 
3D space on the microscopic level without any adjustable parameters. The calculation technique for r\ is 
similar to that proposed in the case of the quartic distribution (n]]. New special functions appearing in 
the construction of the phase-transition theory using the sextic distribution were considered in • In 
the case of the p 6 model, we exploit the special functions with two arguments more complicated than the 
parabolic cylinder functions with one argument for the p 4 model. 



3. Basic relations 



We consider a 3D Ising-like system on a simple cubic lattice with N sites and period c. The Hamilto- 
nian of such a system has the form 



(3.1) 



where si is the operator of the z-component of spin at the ith site, having two eigenvalues +1 and -1. The 
interaction potential is an exponentially decreasing function 



0(rij) = Aexp(-rij/fo). 



(3.2) 



Here, A is a constant, ry is the interparticle distance, and b is the radius of effective interaction. In the 
representation of the CV pk, the partition function of the system in the absence of an external field J€ 
{h - \i^J€ = 0, pb is the Bohr magneton) can be written in the form 



Z = 2"2 



AT 9 (iVi-l)/2 



QolQlP)] 



Ni 



(ex P J-i £ [d> 



(fc)-d'(Bi,B')]p k p_ k 



(l + Ag-i — )exp 

h mv.Ni- 1 ^..x^ 



E Pkr--Pk 2 ,<5k 1+ -+k 2/ 



(dp) 



(3.3) 



where B 1 = B' I s, N\ = JV's" 3 , B' = [b \/2] \ N' = Ns~ 3 , s = BIB', B = nlc is the boundary of Brillouin 
half-zone, and 5k 1+ - +k4 is the Kronecker symbol. For the coefficient d'{k), we have 



d'{k) = a' 2 -pci)(k). 



(3.4) 



Here, ft = \l{kT) is the inverse temperature. For the Fourier transform of the interaction potential, we 
use the following approximation jiillill: 



<l(fc) = 



<l>(0)(l-2/7 2 fc 2 ), k^B', 
0, B'<k^B. 



(3.5) 



The quantities Qo, Q{P) and R 2 i are ultimate functions of the initial coefficients a' 2l {l — 0,1,2,3) 12111611. 
The coefficients a' 2l are determined by special functions of two arguments and are dependent on the 
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ratio of the effective interaction radius b to the lattice constant c, i.e., on the microscopic parameters of 
the system (see, for example, Iil711). 

The correction, which is introduced by the operator A g , is considered in the linear approximation in 

A<£(fc) = q-2b 2 /3®(0)k 2 . (3.6) 

The quantity A<t(k) corresponds to the deviation (5<b(k) from the average value f)<b(Bi,B'). Here, q = 
qp<t(0), q defines the geometric mean value of k 2 on the interval (1/5,1]. In the above-mentioned ap- 
proximation, we arrive at the expression 

-m 1 \4C{h,a)] 2 d 6 . , „ 

4 = r E o —{NT 4 E 

x ^exp[-i(ki+k 2 + k 3 -i-k)li-i(k4 + k5+k 6 -k)l2], (3.7) 
li,l 2 



which defines the operator A g accurate to within the term proportional to g Pk (the terms propor- 

tional to the higher orders of operators d/dp^ are not taken into account). The summation over the sites 
li, I2 in <3.7) is carried out for the lattice with period d = nb\/2. The role of Ag(k) is played by the quantity 

Ag(fc) = , (3.8) 

5 l-S 2 (2^)- 2 A<P(fc) 

where S2 = {2n) 2 [241 a^) 112 F2(h, a). The forms of the functions C{h,a) and F2{h,a) as well as of their ar- 
guments h and a are presented in the next section. We assume that Ag operates only on 
exp(-ii?2ZfcsB 1 PkP-k) m <3.3> . This assumption is associated with small contributions from R4 and 
Re in comparison with the contribution from R2 (in particular, i? 4 / [6R 2 ) ~ 10~ 4 jllll ). 



4. Partition function of the system with allowance for the correction for 
the potential averaging 

The integration over the zeroth, first, second, . . . , nth layers of the CV phase space leads to the rep- 
resentation of the partition function in the form of a product of the partial partition functions Q n of 
individual layers and the integral of the "smoothed" effective distribution of fluctuations. We have 



Z = 2"2 



iV 9 (JV„ +1 -l)/2 



QoQi---QnlQiPn)] N " +l J r 6 (n+1 \p)(dp) N " +1 . (4.1) 



The effective sextic distribution of fluctuations in the {n + l)th block structure is written as follows: 

(4.2) 



# B (n+1) (p) = exp 



3 



■\ Z 4+i(fc)PkP-k-£ ^ t-1 E Pkr--Pk 2/ 5k 1+ - + k 2( 



Here, B n+ i = B' s~ (n+1 \ N n+ i = N's- 3{n+1) , and d n+ i(k), df +l \ a^ +1) are the renormalized values of 
coefficients d'[k), a' v a' 6 after integration over n + 1 layers of the phase space of CV. 

In comparison with the results obtained earlier without taking into account the correction for the 
potential averaging (see, for example, |2j,|8l|l6j,|l7|]), the expression d4.lt includes new quantities. In par- 
ticular, the function 

/ {h n , a n ) = -48s 9l2 C m {h n , anWlinn.Zn) -^=^0 (4.3) 



appearing in Q n characterizes the correction to the partial partition functions. The basic arguments 
h n - d~n {Bn+i:Bn) [61 a^) 112 and a n = \/6a^ /[l5(a^"') 3/2 ] are expressed in terms of the coefficients 
d n {B n+ \,B n ) (the average value of d n {k) in the nth layer of the CV phase space), a^ 1 and flg . The in- 
termediate variables n n = (6s 3 ) 1 ' 2 F 2 {h n ,a n )/C ll2 (h n ,a n ) and£„ = \/6s~ 3l2 N{h n ,a„)/[l5C 3l2 (h n ,a n )] 
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are functions of h n and a n . The special functions C{h n ,a n ) — —F^\}i n ,OL n )+3F|(h„,a„) and N{h n ,a„) = 



nt&n) ~ ]-5F4{hfi,(Xn)F2[hfi t (Xfi) + 30F~ [h n , oc n ) are combinations of the functions F2ii.h n ,(Xn) 



Izi{h n ,a n )/ Io{h n ,a n ), where l2i(h n ,a n ) = J °° t 21 exp [-h n t 2 - t A - a n t 6 ) dt. The relation <4.3> for 
f(h n ,a n ), except 

1 + a' 1 + a' 1 + a' , 

2 ° * 5^, (4.4) 

contains the factor 



, £^^_ (45) 

" " 24 F 2 [h n> a n ) l + a' Q l + a[ l + a'^ /3O(0) ' 

For the quantity a' (1 , which defines the correction for the averaging of the Fourier transform of the poten- 
tial in the nth layer of the phase space of CV, we obtain 

a' n = U4n 2 s 6 F*{T] n ,{, n )qt n S8 Q . (4.6) 

The expressions for 2Pq (appearing in f(h n ,a„)), (appearing in t n ) and 3&q (appearing in a' n ) are 
presented in |ilE3l- At large values of the RG parameter s, there emerge large intervals of wave vectors, 
in which <J>(fc) is averaged. In this case (s > 5), the correction /3<b{k) - /3<b{B n+ \,B n ) is substantial, so that 
its accounting in the linear approximation is incorrect. 

5. Analysis of recurrence relations for the p 6 model. Critical exponent 17 

The coefficients d n+ i(k), a[ n+1) and a l 6 n+1) in {4^ satisfy the following RR: 

dn+i(B n+2 ,B n+ i) = (ay) Y(h n ,a n )- q — 5 2 2 — 1 T~ 



= afs- 3 B(h n ,a n ), 
gtn+l) ' "^ 3/2 

"6 



s- 6 D(^ n ,a„). (5.1) 



The contributions to the functions 

Y(h n ,a n ) = Y(h n , a„) [1 - G(h n , a n )jJ \ , 

B{h n ,a n ) = B(h n ,a n ) [1 +J^(fi„,a„) < g , ]. 

D{h n ,a n ) = D(/j„,a„)[l-if(/i„,a„)@o] (5.2) 

from the potential averaging are given by the terms G(/z n , a n )^o» ^{h n , ocn^a and J£(h n , a n )2>o- Here, 

r(fc„,a„) = s 3/2 F 2 (?? n ,f„)C- 1/2 (/i„,a„), 

B(.h n ,a n ) = s 6 C(.T] n ,S n )C- l (.h n ,a n ), 

D{h n ,a n ) = s 2V2 N{T] n ,( n )C- 3l2 {h n ,a n ), (5.3) 

and 



G{h n ,a n ) = 288s y " L F 2 '(T ln ,Z n )C LU (h n ,a n )- 
X[h n ,a n ) = l728s 3 ' 2 ^4^C m (h n ,a n )- qtn 



X(h n ,a n ) = 172805- 3/2 F ^ n 'l n ] C m (h n ,a n )^. (5.4) 
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The quantities si§, <£q and @o as well as the quantities &q and SSq mentioned above appear due to the 
inclusion of the averaging correction. They are calculated with the help of the summation over the dis- 
tances to the particles located at the lattice sites (see |i|,[ll|]). In terms of the variables 



- — " —d n (.0), 

1 + a' 1 + a\ ! + 



a[ n) 



(l + O^l + a'J 2 (l + a'„_J 2 4 



s 6 s 6 s 6 



w„ = : ~ - 7 - 7 4 n) , (5.5) 



the RR <5.1) assume the forms 



f n +i = r [-q+[u n r Y(h n ,a n )] 

1 + a' L J 



s - 

Wn+i = — ~{u n f l2 b(h„,a n ). (5.6) 



There are two essential distinctions between the RR H5.6) and those obtained without involving the 
correction for A<t>(fc) The first of them consists in the specific substitution of variables <5.5> . 

which differs from the corresponding substitution without the correction by including the factors 
(l + af)) [l + a' 1 )---[l + a' n _ 1 ). The second distinction concerns the transformation of special functions 
Y{h n , a n ), B(h n ,a n ) and D{h 

n,(x>n) (5-3) into the functions Y{h n ,tx n ), B{h nr tx n ) and D{h n , oc n ) iS.2) . This 
distinction is associated with a shift of the fixed-point coordinates and with corrections to the critical 
exponents of thermodynamic functions. 

A particular solution ofRR <5.6t is a new fixed point (r, u, w), which differs at A<J>(fc) + from the fixed 
point ( r (0) , u (0) , u/ (0) ) for the case of AO>(A;) = dEIES. The coordinates of the fixed point of RR {5jD can 
be expressed as follows: 



Here, 



f = -/>t>(0), u = (p[p<S>m, w> = y[)6<P(0)f. (5.7) 



/ = q^Y{h,a)-hl v^] ^Y {h,a)-{\ + a m )hj {s 2 v^)] \ 

<p = q 2 [l-{l + a' (0) )s- 2 ] 2 \Y{h,a)-{l + a m )h/(s 2 V6}]~ 2 , 

iff = {l + a m y\(pf l2 b{h,a), (5.8) 



and 



h = V& 



r+ q 



a 



\/6 w 



15 [uf 12 ' 
cc m = lUn 2 s 6 F 4 2 {q {0 \^)qt m SS , 

f (0) = t n (u {0 \h {0 \a m ). (5.9) 

The quantities h m , a t0) and £ (0) describe, respectively the basic (h n , a„) and intermediate (x\ n , S, n ) 
arguments at the fixed point obtained without involving the correction for the potential averaging. In the 
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linear approximation in A<P(fc), we have 

/ = fJi + { \y' h {h {0 \ a m ) h m /V&-Y {h m ,a {0) ) / Ve] Ah + Y' a (h®, a m ) h m / V6Aa 

- Y [h m , a m ) G [h {0 \ a (0 Vo^°V (l - s~ 2 ) { [ Y {h m , a m )- h m / Vs] 

x [ y [h m ,a m ) - h m /{s 2 V6)]y 1 + a m h™ /{s 2 V&)\y {h [Q \a {0) ) - h m /{s 2 v^)]" 1 ), 

<p = <p \\+2{-\Yl{h {0 \a {0) )-\/{s 2 V^\Ah-Y' a ^ 

+ a m h m /{s 2 V6)\\y {h m , a m ) - h m /{s 2 v^] - 2a m s~ 2 /(l - s" 2 )) , 

-3f{ [ Y (fc (0) , a (0) ) - h m /{s 2 yfi)\ D' h {h m , a (0) ) /3 - (fc (0) , a (0) ) - l/{s 2 Vg)] 

><D{h^\a m )\Ah + {\Y{h^\a m )-h^y{s 2 ^]D' a {h m ,a m )/3 

- Y' a {h w ,a m )D{h {0 \a m )}Aa+Y {h m ,a m )D {h m , a (0) ) G {h m , a'°Vo 

+ a' m h m D{h {0 \a m )/{s 2 V6)){ \Y{h (0 \a m )- h {0) /{s 2 Vg)] D^'V)}" 1 



(5.10) 



The quantities fy, (po and y/o characterize the fixed-point coordinates in the case when r\ = 0. Formulas for 
the derivatives Y' h (fa (0) ,a t0) ), Y^{hP\a m ), D' h {h [a \a m ) and D' a {h m ,a {0} ) can be written using series 
expansions of the corresponding functions in the vicinity of the fixed point (3- The differences Ah = 
h - fo (0) and Aa - a - a (0) determine the displacements of the basic arguments h„ and a n at the fixed 
points (f , u, w) and (r t0) , i/ (0) , w m ). 

The RR i5.6\ make it possible to find the elements of the RG linear transformation matrix. These matrix 
elements Bh (i = 1,2,3; j = 1,2,3) can be presented in the following forms (the linear approximation in 
A<5(fc)): 

Rn = Rn{l-a m ) + RU h) bh + RU a) Aa + R™s*o, 
R 22 = R 22 {l-2a m ) + R 2 1 2 h) Ah + R z l 2 a) Aa + R {2 2 )c £o, 
£33 = R 3 3{l-3a' m ) + R ( 3 \ h) Ah + R ( 3 \ a) Aa + R (2 : j@ , 
Rij = Rf/m li ~ JV2 , 

- jStl-a^+JiSfAh + ligfAa + Ugflb. ^3 = 1,2. (5.11) 

It should be noted that formulas for the quantities Ru and By (i + j) from dS.llt coincide with the 
corresponding expressions for the matrix elements obtained without taking into account the correction 
for the potential averaging |2|. The contributions to the matrix elements Rij from terms R^ Ah and 

J IJ 

Act correspond to a fixed-point shift due to the inclusion of the dependence of the Fourier transform 
of the interaction potential on the wave vector. The terms like i?! 2 '^o» R^^o and R^3>o describe a direct 
contribution to Bu from the correction for averaging. 
The explicit solutions of RR 

f„ = T + ciE^ + C2wf 2 \m~ ll2 E^ + C3wf 3 \uy l E^, 
u n = u + ciw^{u) 1/2 E 1 +c 2 E 2 +c 3 w 2 ° 3 ) {ur 1/2 E", 

w n = w+ciw 3 ®uE 1 +c 2 wf 2 \u) U2 E 2 +c 3 E 3 1 (5.12) 
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in terms of HS.St read 



d n {B n+ \,B n ) - s 



-An) 



-In 



n-1 

II ( 1 + a 'm) 
m=0 



n-1 



n(i+o 2 



m=0 



-6;? 



m=0 



f + q + ci£f + c 2 w^{u)~ ll2 E^ + c 3 (m) -1 ^ 



[fi + Cj U>™ (W) 1/2 £i" + C 2 £ 2 " + C 3 U>™ (")" 1/2 4"] - 

H> + ci h>^ fiEf + c 2 (w) 1/2 £ 2 " + c 3 E$] . (5.13) 



Here, f, ft and w are given in <5.7> and d5.10t . The coefficients q are obtained from the initial condi- 
tions at n — 0. The temperature-independent quantities w.) determine the eigenvectors of the RG linear 
transformation matrix, and £/ are the eigenvalues of this matrix. 

The main distinction between the solutions <5.13> and the solutions of RR in the absence of the cor- 
rection for the potential averaging | 2, 8] lies in the availability of factors of the type (l + a' m ). Taking into 
account the fact that linim^oo a' m {T c ) = a' (0) at T = T c , we obtain the following asymptotics in n for the 
quantities d n , a[ n) and a ( 6 n) from H5.13) : 



d n {Bn+i,B n ) 



{f+q)s 



-n{2-n) 



US 



= ws 



-2n(2-j)) 
-3n(2-)7) 



The quantity r\ is given by the formula 



T) = 



a*™ 
Ins 



(5.14) 



(5.15) 



and corresponds to the critical exponent of the correlation function. Thus, the correction for the potential 
averaging being involved in calculating the partition function of the system, leads to a change of the 
asymptotics for the coefficients d n , a 4 and al at T = T c (in contrast to the case of A<P(fc) = 0, the 
exponents of these coefficients contain the quantity rj). 

The renormalization of the critical exponent of the correlation length v = lns/ln£i, compared to the 
case of rj = 0, is associated with a change of the larger eigenvalue (£i > 1) of the RG linear transformation 
matrix. In contrast to v, the critical exponent of the susceptibility y = (2 - 77) v explicitly depends on rj. The 
specific heat of the system is characterized by the exponent a — 2 - 3v, the expression for which contains 
a renormalized critical exponent of the correlation length v. 



6. Discussion and conclusions 

Within the CV approach, an analytic method for calculating the free energy of a 3D Ising-like system 
near the critical point was elaborated in jllll with the allowance for the simplest non-Gaussian fluctu- 
ation distribution (the p 4 model) and the critical exponent 77. The critical exponent of the correlation 
function rj = 0.024 was found in |11] by including the correction for the potential averaging in the course 
of calculating the partition function of the system. This value of rj is in accord with the other authors' 
data. For comparison, the exponents 77 = 0.0335(25), 77 = 0.0362(8) and 77 = 0.033 were obtained within 
the framework of the field-theory approach (7-loop calculations) 1 19], Monte Carlo simulations 1 20] and 
non-perturbative RG approach (the order <5 4 of the derivative expansion) llill . respectively. Some differ- 
ence between the value of 77 = 0.024 and other authors' data can be connected with the approximations 
in calculations, in particular, with the approximations which are made within the CV method (the sim- 
plest non-Gaussian distribution is used for obtaining 77 = 0.024; the correction inserted by the operator 
Ag is considered in the linear approximation in A<J>(fc); the terms proportional to the higher orders of 
operators d/dpk are not taken into account in the expression for A g ; the operator A g in <3.3> acts only 
on the first term in the exponent). As is established in |11Q , the inclusion of a nonzero exponent 77 within 
the CV method leads to a reduction of the value of the critical exponent for the correlation length, v (in 
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comparison with the case of 77 = 0). For better quantitative estimates of v and other renormalized critical 
exponents, it is necessary to use the non-Gaussian approximations higher than the quartic distribution 
(e.g., the sextic one). 

This paper supplements the previous study 1 11] based on the p 4 model. In the present publication, the 
correction for the averaging of the Fourier transform of the interaction potential is taken into account 
using the sextic fluctuation distribution (the p 6 model). The effect of the mentioned correction on the 
critical behaviour of a 3D uniaxial magnet is investigated in the linear approximation. 

Extending the method for the layer-by-layer integration of the partition function of the system to the 
case of the p 6 model, the RR for the coefficients of the effective distributions are written and analysed. 
It is shown that the inclusion of the correction for the potential averaging gives rise to a change of the 
asymptotics for the RR solutions at T - T c . 

The critical exponents of the correlation length, susceptibility and specific heat are renormalized due 
to the above-mentioned correction. 

Our analytic representations acquire a more general and complete character as compared with the 
case when r\ — 0. An explicit expression (4.1) for the partition function allows us to calculate the free 
energy and other thermodynamic characteristics of the system near T c taking into account the non- 
Gaussian sextic distribution and the small critical exponent r\. 

An analytic procedure for the calculation of the critical exponent of the correlation function devel- 
oped in this paper on the basis of the p 6 model for a one-component spin system may be generalized to 
the case of a system with an rc-component order parameter. 
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KpnTMMHa noBefliHio Tpi/iBi/iMipHOi i3i/iHronoAi6HOi'ci/iCTeMi/i 
b Ha6/in>i<eHHi MOfle/ii p 6 : Po/ib nonpaBKi/i Ha ycepeflHeHhm 
noTemjiasiy 

LB. nn/iK)iP, M.B. y/ifli<P 

iHCTHTyT (|)i3HKW KOHfleHCOBaHWX CUCTeM HAH yKpai'HM, By/1. CBEHL4iL4bKOrO, 1, 7901 1 ZlbBiB, YKpaiHa 

flepxaBHwR Bwinnti HaB4a/ibHwi7i 3aK/iafl 'TipHUHO-eKOHOMiHHUii KO/ieflx", 
By/i. Cryca, 17, 80100 HepBOHOrpafl, YKpaiHa 

3 BMKopucTaHHflM MeTOfly KO/ieKTHBHUx 3MiHHHx (K3) BHB4aeTbCfl KpnTUHHa noBefliHKa cudeM, AKi Ha/iexaTb 
flo K/iacy yHiBepca/ibHOCTi TpuBWMipHOi' MOfle/ii l3WHra. CTaTHCTi/iHHa cyMa OflHOKOMnoHeHTHOi cniHOBOi' cw- 
CTeMM o64nc/iK)eTbCfl lu/ihxom iHTerpyBam-ifl 3a luapaMU cf>a30Boro npodopy K3 b Ha6/inxeHHi HeraycoBoro 
LuecTupHoro po3nofli/iy cp/iyKTyai4iCi napaMeTpa nopaflKy (MOfle/ib p 6 ). Oco6/iwBicTO 3anponoHOBaHoro po3pa- 
xyHKy e npuwHHTTO flo yBam 3a/iexHOCTi (f>yp'e-o6pa3y noiei-mia/iy B3aeMOfliT Bifl XBW/ibOBOro BeKTopa. Bpa- 
xyBaHHa nonpaBKW Ha ycepeflHem-ifl noTei-mia/iy npwBOflHTb flo BiflMiHHOro Bifl Hy/w KpnTWHHOro noKa3HWKa 
Kope/wi4iPiHOi' c(>yHKL4M rj i nepeHopMyBam-ifl 3Ha4eHb imunx KpnTW4Hnx noi<a3HHKiB. Bw,qi/ieHO BHecKH Bifl u,\a 
nonpaBKH b peKypeHTHi cniBBiflHOLueHHfl fl/ifl MOfle/ii p 6 , KoopflUHa™ cf>iKC0BaH0i'T04KH Ta e/ieMeH™ MaTpni^i 
/liHMHOro nepeTBopeHHfl peH0pMa/ii3ai_(ii7iH0i' rpynw. Bnpa3 fl/ia Ma/ioro KpwTH4Horo noKa3HWKa rj OTpuMaHO y 
BumoMy HeraycoBOMy Ha6/iw>KeHHi. 

K/iHJMOBi c/ioBa: rpiABViMipna i3HHronofli6Ha cmcieMa, KpvirvmHa noBefliHKa, wecmpHHfi po3nofli/i, 
ycepeflHeHHH noTeni^iany, Maying KpmmHnti noKd3Hi/iK 
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